Association of the FGF4L2 retrogene with fibrocartilaginous embolic myelopathy in dogs.

BACKGROUND: Fibrocartilaginous embolic myelopathy (FCE) is a well-documented condition in dogs although rarely reported in chondrodystrophic breeds. Genetic associations have not been defined. OBJECTIVES: Define the association of the chondrodystrophy-associated FGF4L2 retrogene with histopathologically confirmed cases of FCE. ANIMALS: Ninety-eight dogs with a histopathologic diagnosis of FCE. METHODS: Retrospective multicenter study. Dogs were genotyped for the FGF4L2 and FGF4L1 retrogenes using DNA extracted from formalin-fixed, paraffin-embedded tissue. Associations between breed, FCE and retrogene status were investigated with reference to a hospital population and known breed and general population allele frequencies. RESULTS: FGF4L2 genotype was defined in 89 FCE cases. Fibrocartilaginous embolic myelopathy was present in 22 dogs from FGF4L2-segregating breeds with allele frequencies of ≥5%; however, all dogs were wild type. Two Labrador retrievers with FCE carried FGF4L2 alleles. Frequency of the FGF4L2 allele was significantly (P < .001) and negatively associated with FCE relative to predicted hospital-population dogs. FCE was overrepresented in Boxer, Great Dane, Yorkshire Terrier, Bernese Mountain Dog, Miniature Schnauzer, Rottweiler, and Shetland Sheepdog breeds. CONCLUSIONS AND CLINICAL IMPORTANCE: Study data based on genotypically and histopathologically defined cases support the historical observation that FCE is uncommon in chondrodystrophic dog breeds. FGF4 plays an important role in angiogenesis and vascular integrity; anatomical studies comparing chondrodystrophic and non-chondrodystrophic dogs might provide insight into the pathogenesis of FCE.

Naturally occurring canine laminopathy leading to a dilated and fibrosing cardiomyopathy in the Nova Scotia Duck Tolling Retriever.

Dilated cardiomyopathy (DCM) is characterized by decreased systolic function and dilation of one or both ventricles, often leading to heart failure or sudden death. Two 10-month-old sibling Nova Scotia Duck Tolling Retrievers (NSDTR) died acutely with evidence of dilated cardiomyopathy with myocardial fibrosis. Association analysis using two cases and 35 controls identified three candidate regions homozygous in the two cases. Whole genome sequencing identified a frameshift deletion in the LMNA gene (NC_049228.1:g.41688530del, NP_001274080:p.(Asp576ThrfsTer124)). Three retrospectively identified NSDTRs with sudden death before 2 years of age and severe myocardial fibrosis were also homozygous for the deletion. One 5 year old with sudden death and myocardial fibrosis was heterozygous for the deletion. This variant was not identified in 722 dogs of other breeds, nor was it identified to be homozygous in 784 NSDTR. LMNA codes for lamin A/C proteins, which are type V intermediate filaments that provide structural support to the nuclear membrane. In humans, LMNA variants can cause DCM with sudden death as well as diseases of striated muscles, lipodystrophy, neuropathies, and accelerated aging disorders. This frameshift deletion is predicted to affect processing of prelamin A into lamin A. Pedigree analysis in the NSDTR and functional evaluation of heterozygotes is consistent with a predominantly recessive mode of inheritance and possibly low penetrance in heterozygotes in contrast to people, where most pathogenic LMNA variants are dominantly inherited.

A variant in the 5'UTR of ERBB4 is associated with lifespan in Golden Retrievers.

Genome-wide association studies (GWAS) in long-lived human populations have led to identification of variants associated with Alzheimer's disease and cardiovascular disease, the latter being the most common cause of mortality in people worldwide. In contrast, naturally occurring cancer represents the leading cause of death in pet dogs, and specific breeds like the Golden Retriever (GR) carry up to a 65% cancer-related death rate. We hypothesized that GWAS of long-lived GRs might lead to the identification of genetic variants capable of modifying longevity within this cancer-predisposed breed. A GWAS was performed comparing GR dogs ≥ 14 years to dogs dying prior to age 12 which revealed a significant association to ERBB4, the only member of the epidermal growth factor receptor family capable of serving as both a tumor suppressor gene and an oncogene. No coding variants were identified, however, distinct haplotypes in the 5'UTR were associated with reduced lifespan in two separate populations of GR dogs. When all GR dogs were analyzed together (n = 304), the presence of haplotype 3 was associated with shorter survival (11.8 years vs. 12.8 years, p = 0.024). GRs homozygous for haplotype 3 had the shortest survival, and GRs homozygous for haplotype 1 had the longest survival (11.6 years vs. 13.5 years, p = 0.0008). Sub-analyses revealed that the difference in lifespan for GRs carrying at least 1 copy of haplotype 3 was specific to female dogs (p = 0.009), whereas survival remained significantly different in both male and female GRs homozygous for haplotype 1 or haplotype 3 (p = 0.026 and p = 0.009, respectively). Taken together, these findings implicate a potential role for ERBB4 in GR longevity and provide evidence that within-breed canine lifespan studies could serve as a mechanism to identify favorable or disease-modifying variants important to the axis of aging and cancer.

Ancient segmentally duplicated LCORL retrocopies in equids.

LINE-1 is an active transposable element encoding proteins capable of inserting host gene retrocopies, resulting in retro-copy number variants (retroCNVs) between individuals. Here, we performed retroCNV discovery using 86 equids and identified 437 retrocopy insertions. Only 5 retroCNVs were shared between horses and other equids, indicating that the majority of retroCNVs inserted after the species diverged. A large number (17-35 copies) of segmentally duplicated Ligand Dependent Nuclear Receptor Corepressor Like (LCORL) retrocopies were present in all equids but absent from other extant perissodactyls. The majority of LCORL transcripts in horses and donkeys originate from the retrocopies. The initial LCORL retrotransposition occurred 18 million years ago (17-19 95% CI), which is coincident with the increase in body size, reduction in digit number, and changes in dentition that characterized equid evolution. Evolutionary conservation of the LCORL retrocopy segmental amplification in the Equidae family, high expression levels and the ancient timeline for LCORL retrotransposition support a functional role for this structural variant.

FGF4L2 retrogene copy number is associated with intervertebral disc calcification and vertebral geometry in Nova Scotia Duck Tolling Retrievers.

OBJECTIVES: To evaluate the effects of the chondrodystrophy-associated FGF4L2 retrogene on intervertebral disc (IVD) calcification and vertebral geometry. ANIMALS: 22 Nova Scotia Duck Tolling Retrievers (NSDTR) with no FGF4L2 retrogene (n = 7, wild-type dogs), 1 retrogene copy (8, heterozygous dogs), or 2 retrogene copies (7, homozygous dogs). PROCEDURES: Computed tomography (CT) scans of the vertebral column were analyzed using computer-aided design (CAD) software. IVD calcification, vertebral column length, and vertebral geometry of the third cervical (C3), 13th thoracic (T13), and first lumbar (L1) vertebrae were compared. RESULTS: IVD calcification was not found in wild-type dogs. IVD calcification was more frequent in homozygous dogs than heterozygous (P = .008) or wild-type dogs (P < .001) and in heterozygous dogs compared to wild-type dogs (P < .001). Four IVDs were subclinically herniated in 3 dogs (2 homozygous, 1 heterozygous). Calcified IVD had a greater volume and surface area in heterozygous dogs than homozygous dogs. C3 vertebral canal height-to-width ratio was greater in homozygous dogs than heterozygous dogs (P = .044) and wild-type dogs (P = .010). CLINICAL RELEVANCE: IVD calcification and vertebral geometry can be analyzed using CAD software. The presence of 1 or 2 FGF4L2 copies in the absence of the FGF4L1 retrogene has an additive effect on the number of calcified IVD and a minor effect on vertebral geometry in NSDTR dogs. Data support the use of FGF4L2 phenotyping to reduce clinical disease in segregating breeds and to monitor the introduction of wild-type alleles into fixed breed populations.

Identification of Genetic Risk Factors for Monogenic and Complex Canine Diseases.

Advances in DNA sequencing and other technologies have greatly facilitated the identification of genetic risk factors for inherited diseases in dogs. We review recent technological developments based on selected examples from canine disease genetics. The identification of disease-causing variants in dogs with monogenic diseases may become a widely employed diagnostic approach in clinical veterinary medicine in the not-too-distant future. Diseases with complex modes of inheritance continue to pose challenges to researchers but have also become much more tangible than in the past. In addition to strategies for identifying genetic risk factors, we provide some thoughts on the interpretation of sequence variants that are largely inspired by developments in human clinical genetics.

Ocular morphologic traits in the American Cocker Spaniel may confer primary angle closure glaucoma susceptibility.

Acute primary angle closure glaucoma is a potentially blinding ophthalmic emergency requiring prompt treatment to lower the elevated intraocular pressure in humans and dogs. The PACG in most of canine breeds is epidemiologically similar to humans with older and female patients overrepresented with the condition. The American Cocker Spaniel (ACS) is among the most common breeds observed with PACG development in dogs. This study initially sought to identify genetic risk factors to explain the high prevalence of PACG in ACSs by using a case-control breed-matched genome-wide association study. However, the GWAS failed to identify candidate loci associated with PACG in this breed. This study then assessed intrinsic ocular morphologic traits that may relate to PACG susceptibility in this breed. Normal ACSs without glaucoma have a crowded anterior ocular segment and narrow iridocorneal angle and ciliary cleft, which is consistent with anatomical risk factors identified in humans. The ACSs showed unique features consisting of posterior bowing of iris and longer iridolenticular contact, which mirrors reverse pupillary block and pigment dispersion syndrome in humans. The ACS could hold potential to serve as an animal model of naturally occurring PACG in humans.

A genome-wide association study to investigate genetic loci associated with primary glaucoma in American Cocker Spaniels.

OBJECTIVE: To identify genetic associations with primary glaucoma (PG) in American Cocker Spaniels using a genome-wide association study (GWAS). ANIMALS: A nationwide ambidirectional case-control cohort study was performed in American Cocker Spaniels that had an ophthalmic examination performed by a veterinarian. Ninety-four dogs with PG (cases) and 111 dogs without glaucoma (controls) met phenotypic criteria and had a blood sample collected after receiving informed owner consent. PROCEDURES: Genomic DNA was extracted from whole blood samples and genotyped (CanineHD BeadChip, Illumina Inc). A case-control GWAS using a linear mixed model was performed, and 3 significance thresholds were calculated (1) using a Bonferroni correction on all single nucleotide polymorphisms (SNPs) included in the GWAS, (2) using a Bonferroni correction on only the unlinked SNPs from a pruned data set, and (3) using 10,000 random phenotype permutations. RESULTS: Following genotype data quality control, 89 cases and 93 controls were included in the GWAS. We identified an association on canine chromosome (CFA10); however, it did not reach statistical significance. Potential candidate genes within the surrounding linkage disequilibrium interval include coiled-coil domain containing 85A (CCDC85A) and extracellular growth factor containing fibulin extracellular matrix protein 1 (EFEMP1). CLINICAL RELEVANCE: Primary glaucoma in the American Cocker Spaniel is a complex heterogeneous disease that may be influenced by a locus on CFA10. The candidate genes CCDC85A and EFEMP1 within the identified linkage disequilibrium interval have been shown to be involved in human open-angle glaucoma.

An SNN retrocopy insertion upstream of GPR22 is associated with dark red coat color in Poodles.

Pigment production and distribution is controlled through multiple genes, resulting in a wide range of coat color phenotypes in dogs. Dogs that produce only the pheomelanin pigment vary in intensity from white to deep red. The Poodle breed has a wide range of officially recognized coat colors, including the pheomelanin-based white, cream, apricot, and red coat colors, which are not fully explained by the previously identified genetic variants involved in pigment intensity. Here, a genome-wide association study for pheomelanin intensity was performed in Poodles which identified an association on canine chromosome 18. Whole-genome sequencing data revealed an SNN retrocopy insertion (SNNL1) in apricot and red Poodles within the associated region on chromosome 18. While equal numbers of melanocytes were observed in all Poodle skin hair bulbs, higher melanin content was observed in the darker Poodles. Several genes involved in melanogenesis were also identified as highly overexpressed in red Poodle skin. The most differentially expressed gene however was GPR22, which was highly expressed in red Poodle skin while unexpressed in white Poodle skin (log2 fold change in expression 6.1, P < 0.001). GPR22 is an orphan G-protein-coupled receptor normally expressed exclusively in the brain and heart. The SNNL1 retrocopy inserted 2.8 kb upstream of GPR22 and is likely disrupting regulation of the gene, resulting in atypical expression in the skin. Thus, we identify the SNNL1 insertion as a candidate variant for the CFA18 pheomelanin intensity locus in red Poodles.

Recent, full-length gene retrocopies are common in canids.

Gene retrocopies arise from the reverse transcription and insertion into the genome of processed mRNA transcripts. Although many retrocopies have acquired mutations that render them functionally inactive, most mammals retain active LINE-1 sequences capable of producing new retrocopies. New retrocopies, referred to as retro copy number variants (retroCNVs), may not be identified by standard variant calling techniques in high-throughput sequencing data. Although multiple functional FGF4 retroCNVs have been associated with skeletal dysplasias in dogs, the full landscape of canid retroCNVs has not been characterized. Here, retroCNV discovery was performed on a whole-genome sequencing data set of 293 canids from 76 breeds. We identified retroCNV parent genes via the presence of mRNA-specific 30-mers, and then identified retroCNV insertion sites through discordant read analysis. In total, we resolved insertion sites for 1911 retroCNVs from 1179 parent genes, 1236 of which appeared identical to their parent genes. Dogs had on average 54.1 total retroCNVs and 1.4 private retroCNVs. We found evidence of expression in testes for 12% (14/113) of the retroCNVs identified in six Golden Retrievers, including four chimeric transcripts, and 97 retroCNVs also had significantly elevated F ST across dog breeds, possibly indicating selection. We applied our approach to a subset of human genomes and detected an average of 4.2 retroCNVs per sample, highlighting a 13-fold relative increase of retroCNV frequency in dogs. Particularly in canids, retroCNVs are a largely unexplored source of genetic variation that can contribute to genome plasticity and that should be considered when investigating traits and diseases.

SLC25A12 Missense Variant in Nova Scotia Duck Tolling Retrievers Affected by Cerebellar Degeneration-Myositis Complex (CDMC).

We investigated two litters of distantly related Nova Scotia Duck Tolling Retrievers (NSDTR), of which four puppies were affected by cerebellar signs with or without neuromuscular weakness. The phenotype was termed cerebellar degeneration­myositis complex (CDMC). We suspected a heritable condition and initiated a genetic analysis. The genome of one affected dog was sequenced and compared to 565 control genomes. This search yielded a private protein-changing SLC25A12 variant in the affected dog, XM_038584842.1:c.1337C>T, predicted to result in the amino acid change XP_038440770.1:(p.Pro446Leu). The genotypes at the variant co-segregated with the phenotype as expected for a monogenic autosomal recessive mode of inheritance in both litters. Genotyping of 533 additional NSDTR revealed variant allele frequencies of 3.6% and 1.3% in a European and a North American cohort, respectively. The available clinical and biochemical data, together with current knowledge about SLC25A12 variants and their functional impact in humans, mice, and dogs, suggest the p.Pro446Leu variant is a candidate causative defect for the observed phenotype in the affected dogs.

The Effects of FGF4 Retrogenes on Canine Morphology.

Two FGF4 retrogenes (FGF4L1 on chromosome 18 and FGF4L2 on chromosome 12) have been identified to cause dwarfism across many dog breeds. Some breeds are nearly homozygous for both retrogenes (e.g., Dachshunds) and others are homozygous for just one (e.g., Beagles and Scottish Terriers). Since most breeds do not segregate both of these retrogenes, it is challenging to evaluate their individual effects on long bone length and body size. We identified two dog breeds selected for hunting ability, the Alpine Dachsbracke and the Schweizer Niederlaufhund, that segregate both of these retrogenes. Using individual measurements of height at the shoulder, back length, head width, thorax depth and width, and thoracic limb measurements, we evaluated the combined effects of FGF4 retrogenes within these breeds. We applied multivariable linear regression analysis to determine the effects of retrogene copy numbers on the measurements. Copy numbers of both retrogenes had significant effects reducing height at the shoulders and antebrachial length, with FGF4L1 having a much greater effect than FGF4L2. FGF4L1 alone influenced the degree of carpal valgus and FGF4L2 alone increased head width. Neither retrogene had an effect on thorax width or depth. Selectively breeding dogs with FGF4L1 and without FGF4L2 would likely lead to a reduction in the FGF4L2-related risk of intervertebral disc herniation while maintaining the reduction in leg length resulting from FGF4L1.

The effect of inbreeding, body size and morphology on health in dog breeds.

BACKGROUND: Dog breeds are known for their distinctive body shape, size, coat color, head type and behaviors, features that are relatively similar across members of a breed. Unfortunately, dog breeds are also characterized by distinct predispositions to disease. We explored the relationships between inbreeding, morphology and health using genotype based inbreeding estimates, body weight and insurance data for morbidity. RESULTS: The average inbreeding based on genotype across 227 breeds was Fadj = 0.249 (95% CI 0.235-0.263). There were significant differences in morbidity between breeds with low and high inbreeding (H = 16.49, P = 0.0004). There was also a significant difference in morbidity between brachycephalic breeds and non-brachycephalic breeds (P = 0.0048) and between functionally distinct groups of breeds (H = 14.95 P < 0.0001). Morbidity was modeled using robust regression analysis and both body weight (P < 0.0001) and inbreeding (P = 0.013) were significant (r2 = 0.77). Smaller less inbred breeds were healthier than larger more inbred breeds. CONCLUSIONS: In this study, body size and inbreeding along with deleterious morphologies contributed to increases in necessary health care in dogs.

Dog colour patterns explained by modular promoters of ancient canid origin.

Distinctive colour patterns in dogs are an integral component of canine diversity. Colour pattern differences are thought to have arisen from mutation and artificial selection during and after domestication from wolves but important gaps remain in understanding how these patterns evolved and are genetically controlled. In other mammals, variation at the ASIP gene controls both the temporal and spatial distribution of yellow and black pigments. Here, we identify independent regulatory modules for ventral and hair cycle ASIP expression, and we characterize their action and evolutionary origin. Structural variants define multiple alleles for each regulatory module and are combined in different ways to explain five distinctive dog colour patterns. Phylogenetic analysis reveals that the haplotype combination for one of these patterns is shared with Arctic white wolves and that its hair cycle-specific module probably originated from an extinct canid that diverged from grey wolves more than 2 million years ago. Natural selection for a lighter coat during the Pleistocene provided the genetic framework for widespread colour variation in dogs and wolves.

Correction: Weich, K., et al. Pigment Intensity in Dogs Is Associated with a Copy Number Variant Upstream of KITLG. Genes 2020, 11, 75.

The authors wish to make the following corrections to this paper [...].

Variants at the ASIP locus contribute to coat color darkening in Nellore cattle.

BACKGROUND: Nellore cattle (Bos indicus) are well-known for their adaptation to warm and humid environments. Hair length and coat color may impact heat tolerance. The Nellore breed has been strongly selected for white coat, but bulls generally exhibit darker hair ranging from light grey to black on the head, neck, hump, and knees. Given the potential contribution of coat color variation to the adaptation of cattle populations to tropical and sub-tropical environments, our aim was to map positional and functional candidate genetic variants associated with darkness of hair coat (DHC) in Nellore bulls. RESULTS: We performed a genome-wide association study (GWAS) for DHC using data from 432 Nellore bulls that were genotyped for more than 777 k single nucleotide polymorphism (SNP) markers. A single major association signal was detected in the vicinity of the agouti signaling protein gene (ASIP). The analysis of whole-genome sequence (WGS) data from 21 bulls revealed functional variants that are associated with DHC, including a structural rearrangement involving ASIP (ASIP-SV1). We further characterized this structural variant using Oxford Nanopore sequencing data from 13 Australian Brahman heifers, which share ancestry with Nellore cattle; we found that this variant originates from a 1155-bp deletion followed by an insertion of a transposable element of more than 150 bp that may impact the recruitment of ASIP non-coding exons. CONCLUSIONS: Our results indicate that the variant ASIP sequence causes darker coat pigmentation on specific parts of the body, most likely through a decreased expression of ASIP and consequently an increased production of eumelanin.

Canine DVL2 variant contributes to brachycephalic phenotype and caudal vertebral anomalies.

A frameshift deletion variant in the Wnt pathway gene dishevelled 2 (DVL2) is associated with a truncated, kinked tail ("screw tail") in English Bulldogs, French Bulldogs and Boston Terriers. These breeds are also characterized by distinctive morphological traits, including a wide head, flat face and short-limbed dwarfism, which are characteristic of Robinow syndrome in humans, caused by defects in genes such as DVL1 and DVL3. Based on these phenotypic and genetic similarities, it has previously been hypothesized that the canine DVL2 variant results in a syndromic phenotype called the Robinow-like syndrome. In our study, we investigated the distribution of the DVL2 variant in 1954 dogs from 15 breeds, identifying breeds with allele variation and enabling the dissection of the genotype-phenotype correlation for the first time. With CT examinations in American Staffordshire Terriers, we confirmed that the DVL2 allele is associated with caudal vertebral malformations and a brachycephalic phenotype. We also hypothesize that the variant may be linked to additional health conditions, including brachycephalic obstructive airway syndrome and congenital heart defects. Altogether, our study strengthens the role of DVL2 as one of the contributors to the "bulldog type" morphology and features on the spectrum of human Robinow syndrome.

Special Issue "Molecular Basis of Inherited Diseases in Companion Animals".

The study of inherited diseases in companion animals has exploded over the past 15 years since the publication of the first dog genome in 2005 [...].

A Missense Variant in ALDH5A1 Associated with Canine Succinic Semialdehyde Dehydrogenase Deficiency (SSADHD) in the Saluki Dog.

Dogs provide highly valuable models of human disease due to the similarity in phenotype presentation and the ease of genetic analysis. Seven Saluki puppies were investigated for neurological abnormalities including seizures and altered behavior. Magnetic resonance imaging showed a diffuse, marked reduction in cerebral cortical thickness, and symmetrical T2 hyperintensity in specific brain regions. Cerebral cortical atrophy with vacuolation (status spongiosus) was noted on necropsy. Genome-wide association study of 7 affected and 28 normal Salukis revealed a genome-wide significantly associated region on CFA 35. Whole-genome sequencing of three confirmed cases from three different litters revealed a homozygous missense variant within the aldehyde dehydrogenase 5 family member A1 (ALDH5A1) gene (XM_014110599.2: c.866G>A; XP_013966074.2: p.(Gly288Asp). ALDH5A1 encodes a succinic semialdehyde dehydrogenase (SSADH) enzyme critical in the gamma-aminobutyric acid neurotransmitter (GABA) metabolic pathway. Metabolic screening of affected dogs showed markedly elevated gamma-hydroxybutyric acid in serum, cerebrospinal fluid (CSF) and brain, and elevated succinate semialdehyde in urine, CSF and brain. SSADH activity in the brain of affected dogs was low. Affected Saluki dogs had striking similarities to SSADH deficiency in humans although hydroxybutyric aciduria was absent in affected dogs. ALDH5A1-related SSADH deficiency in Salukis provides a unique translational large animal model for the development of novel therapeutic strategies.